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1.  INTRODUCTION 
4/ 

There  has  been  much  discussion  of  Walsh  spectral  analysis  for  dyadic 

^ - “"-N. 

stationary  processes.  Morettin'L('T97#t  investigated  some  asymptotic  prop¬ 
erties  of  Jthe  finite  Walsh  transforms  of  dyadic  stationary  processes. 

Nagai  %T977T  gave  the  spectral  representations  for  dyadic  stationary  pro¬ 
cesses.  If  we  consider  finite,  dyadic  linear  models  then  the  greatest 
differences  between  dyadic  stationary  processes  and  ordinary  stationary 
ones  appear.  Nagai ‘(T98fT)  and  Nagai  and  Taniguchi  '>(19&4  established  that 
a  dyadic  autoregressive  and  moving  average  o(-MRMAr)  process  of  finite  order 
can  be  expressed  as  a  dyadic  autoregressive  {-BAR)- process  of  finite  order 
and  also  as  a  dyadic  moving  average  {-DMA}  process  of  finite  order.  Nagai 
and  Taniguchi "'(1987)  discussed  the  principal  component  analysis  of  a 

multiple  dyadic  process,  and  also  the  canonical  correlation  analysis. 

- 

Morgettin  4^98+)  gave  a  convenient  survey  for  Walsh  spectral  analysis. 

— Tw  this  papers  consider^  multiple  dyadic  stationary  process  with  ^ 
-the  Walsh  spectral  density  matrix  f  (a),  where  e  is  an  unknown  parameter 
vector.  We  define  a  quasi-maxinium\likelihood  estimator  e  of  e,  and  give 
the  asymptotic  distribution  of  0  under  appropriate  condi tioiis.  In  Section 
3  we  propose  an  information  cri teri on  which  determines  the  order  of  the 
model,  and  show  that  this  criterion  gives  the  consistent  order  estimate. 

In  Section  4,  for  a  finite  order  dyadic  autoregressive  model,  we  propose 
a  simpler  order  determination  criterion,  and  show  that  the  estimated  orderar 
has  strong  consistency.  Also  some  interesting  examples  are  given  in  the 
identification  problem  for  Walsh  spectrum.  In  Section  5  we  consider  a  n 
signal  detection  model  for  dyadic  process  of  finite  order,  and  show  that 

n  / 

this  model  is  equivalent  to  a  dyadic  moving  average  model.  In  Section  6 


. 

;■  O  S  % , 
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we  discuss  a  testing  problem  whether  the  unknown  parameter  e  satisfies 
a  linear  restriction.  Then  we  give  the  asymptotic  distribution  of  the 
likelihood  ratio  criterion  under  the  null  hypothesis. 


J- 


2.  DYADIC  STATIONARY  PROCESSES  AND  ESTIMATION  THEORY 


First  we  introduce  some  basic  concepts  and  notations.  Denote  by  T 
the  set  of  all  nonnegative  integers.  Let  x  and  y  be  nonnegative  real 
numbers  and  have  the  following  binary  expansions: 

oo 

x  =  l  x  21  with  x  =  0  or  1, 

. _  X  X 

00 

y  =  l  y£2*  with  y&  =  0  or  1. 

£=-“ 

Then  the  dyadic  addition^  is  defined  by 

x0y  =  l  |x£-yj2£. 

£=-“ 

A  stochastic  process  (possibly  vector  process)  {Y ( t ) :  t  e  T}  is  said  to  be 
dyadic  stationary  if  the  joint  distribution  of  Y(t-j),  Y ( t^ ) ,  ...»  Y(tn)  is 
the  same  as  that  of  Y ( t-|  ®  t) ,  Y(t2  0  t) ,  . . . ,  Y(tp  0  t)  for  every  finite 
set  of  integers  {t^  tn>  and  for  every  integer  t.  We  denote  by 

(W(t,x) :  0  <_  X  <_  1} ,  t  =  0,1 .  the  system  of  Walsh  functions .  The  prop¬ 
erties  of  Walsh  functions  are  well  known: 

(i)  for  each  t  e  T  and  X  e  [0,1],  the  value  of  W(t,x)  is  only  +1  or  -1, 

(ii )  for  any  s ,  t  e  T, 

W(t,x)W(s,x)  =  W(t0s,  x),  a.e.  x, 

(iii)  for  each  t  6  T  and  x  e  [0,1], 

W(t,x)W(t,y)  =  W(t,  X  >+)  y)  ,  a.e;  y. 

(See  Morettin  (1974).) 

Let  Y(t)  =  ( Y^ (t) ,  ...,  Y  (t)) 1 ;  t  e  T  be  a  q-dimensional  dyadic  stationary 
process  with  zero  mean  and  k-th  order  cumulants 


c 


,  (t, .....  t.  ,)  =  cum{  Y  (t,0tj,  Y  (t90tj,  ..•)  Y  (t.)}, 


)  in  v;  'Wf&inwj.  'irv  ^m'wy.v.'w.v 


t-j ,  •••»t|(_i  6  T.  We  denote  the  covariance  matrices 

r(t,)  =  {c  .  (t.)},  qxq  matrices, 
i  a-|a2  I 

The  statistic 

(H)  N“1 

d(N)(x)  =  l  Y(t)W(t.x) 
t=0 


(2.1) 


is  called  the  finite  Walsh  transform  of  (Y(t):  t  =  0,1,.. .  ,N-1).  Throughout 
this  paper  we  assume  that  N  =  2m,  with  m  a  nonnegative  integer  and  denote 
d^  (x)  =(djN^(x),...  ,d^  (x))  ' .  Here  we  assume  the  fol  lowing. 

ASSUMPTION  1.  For  every  k  and  j  =  1,2 . k-1, 


(2.2) 


tr°  tk_r°  1“  k 


for  all  a^ ,  . . . ,  a^. 

Then  the  Walsh  spectral  density  matrix  and  the  Walsh  cumulant  spectrum 
of  order  k  of  (Y(t)l  are  defined  by 


f(x)  =  l  r(t)W(t,x) , 

t=0 


and 


far..ak(V  xk-l* 


=  I  •••  I  a  (t-i,...,t.  -i)  n  W(t.,X.)  , 

tk_,  1*”  k  1  k  1  j=1  J  J 


k-1 


(2.3) 


respectively.  The  following  proposition  is  due  to  Morettin  (1974), 


PROPOSITION  1.  Under  Assumption  1, 


cum{d^N^ (a-j ) ,  ... ,  d^(Ak)} 
!  k 


,-1 


=  Dn(a1  ©  ...  ©  Ak){fai...ak(Al . Xk-1^  + 


(2.4] 


where  D^(x)  =  £b~Jw(t,x) ,  and  the  term  0(N”^)  is  uniform  with  respect  to 


A1 5  *  •  •  *  Alf.  □ 


Although  we  do  not  assume  the  Gaussianity  of  { Y ( t ) } ,  we  can  make  the 
Gaussian  likelihood  function  L  of  {Y(0),  ....  Y(N-l)},  formally,  and  approxi¬ 
mate  L.  That  is,  we  get 


log  L 


N 

O 

u 


-  1 

{logdet  f_(x)  +  tr  IM(x)f  (a)  }dx  +  constant. 

q  o  n  o 


(2.5: 


where  the  fitted  Walsh  spectral  density  matrix  of  { Y ( t ) )  is  parameterized 

as  f  (a)  ,  e  =  (e. ,  ... ,  e  J  '  6  9=/,  and 
0  ~  I  r 

INU)  =  Fn(x)Fn(x)'  =  Uab(x)},  say, 

1  N-l 

FN(A)  =  -  l  Y{t)W(t,A). 

/N  t=0 

Thus  we  estimate  e  by  the  value  9  which  minimize 


D(f§,IN) 


^{log  det  fg(x)  +  tr  IN(x)fg(x)_1 }dx , 


(2.6 


with  respect  to  9.  Henceforth  we  call  9  quasi-maximum  likelihood  estimator 
of  9.  To  discuss  the  asymptotic  properties  of  9,  the  following  lemma  is  a 
keystone. 


LEMMA  1.  Let  <t>.j(X)  =  {<j>^  (X)},  j  =  l,...,r,  be  qxq  matrix-valued 

Continuous  functions  on  [0,1]  such  that  4> . ( X )  =  4>.(x)'.  Under  Assumption  1 

J  J 

we  can  show  that 

1  rl 

0)  lim  tr  IN(x)<j>  .(x)dX  =  tr  f(X)«t>  .(X) dX ,  in  p,  (2. 

N-^co  Jo  J  JO  J 

(2)  the  quantities 

Aj  =  ^  ftr{(IN(x)  “  f^))+j(^)>dXf  j  =  1 . r, 

have,  asymptotical ly,  a  normal  distribution  with  zero  mean  vector  and 
covariance  matrix  V  whose  (j,m)  element  is 

2  [  tr{f(xU(x)f(xH.(x)}dx 
J  0  J 

q  f1 

+  a  Jr  H-l  H  ^i)(Xl^C){v]f,bc6{X'X^dXd^  <2- 

a  >D ,C ,Q- I  q 


Proof.  Notice  that 

Aj  =  ^  l  f{Iab(x)  '  fab(x),*bi)(>)dx' 
a,b=l  •'0 

By  Proposition  1,  we  have 

E(Iab(x)  -  fab(x))  =  0(1), 

where  0(i)  is  uniform  with  respect  to  X.  Hence  we  obtain 
N 


E(A.)  =  0(N"1/2). 


cum{ I ab ( A ) ,  Icd(y)}  =  ^{cum(da(x),  dd(y ))  cum(db(  a)  ,  dc(y)) 

+  cum(d  (a),  dc(y))  cum(db(x) ,  drf(y)) 
+  cum(da(A),  db(x),  dc(y),  dd(y))} 


1  ,n2 


{D^(X@y)[fad(x)fbc(x)  +  fac(x)fbd(x)] 

+  DN(x©X(±Jy®y)fabcd(X.X,y)  +  0(1)} 


D^(x©X®y®y)  =  D^(0)  =  N, 


we  have 


cum(Aj,Am) 


NfV  l  ^Ja)(x)0^)(y)cun(lab(x),  Icd(y))dX1dX2 

J0J0  a  ,b,c,d=l 

l  f  |’lVJa)(l)*dc)(>')fabcd(i'i"l)dldl) 

a,b,c,d=l  J0-l0 

+  2j1*ba)(>)fac(x)fbd(>)d>[Mo^c)(l')DN(j®'j)dl] 


*  0(r 


DW(X  ©  X)  = 


N  if  | X  ©  y  |  <  1  , 
0  otherwise, 


we  get  that  by  the  continuity  of  4>dc  (m) 


(2.10) 


sC*dc)(“>DN(A®'')d"  ‘ 


^(X)  +  o(l),  if  21  '  X  <  1 


0(1),  otherwise. 


--  »v 


Substituting  the  above  into  (2.10),  we  get 


cum(A.,Am)  =  2 


fl 


tr{f(A)*m(A)f(A)*j(x)}dA 


q 

+  I 

a,b,c,d=1 


1 


4a)(x)<t>dc)(lj)fabcd(x>x’p)dxd,J  +  o(1)-  (2J1) 


Thus  (2.9)  and  (2.11)  imply  our  result  (1).  Also  (2.11)  gives  the  asymptotic 
variance  (2.8).  As  for  the  asymptotic  normality  of  A.,  we  have  only  to 

J 

evaluate  J(J>_3)th  order  cumulant  cum{A.  ,  A.  ,  . . . ,  A .  }.  Here,  without 

1  1 2  J 

loss  of  generality,  we  evaluate  it  for  scalar  process. 

By  Theorem  2.3.2,  p . 21 ,  Brillinger  (1975),  we  have 


cum  1 1  ( A  i )  di  £  ( A  i ) ,  ..., 

=  l  cum(dj.  ( A  j ) ,  (jj)  e  v] )  ...  cum(dj .  {a.j)  ,  (j,i)  e  v$  )  (2.12) 


where  the  summation  runs  over  all  indecomposable  partitions  v  =  VjU  ...  Uvs 
of  the  set  T ( . j ,  i ) ,  j  =  1,2,...,J,  i  =  1,2),  (the  definition  of  indecomposa- 
bility  can  be  found  on  p.20,  Brillinger  (1975)).  By  indecomposabi 1 i ty  of 
the  partitions,  each  vn  contains  at  least  two  elements,  we  have 

S  <  J/2. 

By  Proposition  1,  we  have 


cum(dji  ( Aj) ,  ( j  ,i )  e  v,  )...  cum(d^.(Aj),  (j,i) 


Since 


we  have,  for  J  >  2 


dnU)  = 


if  0  £  A  < 
otherwise, 


i'-f 

Jo 

and  for  J  =  1 


r1  f1 

=  j  ...  j  DN(p1®y2)DN(p20y3)  ...  DN(ys®p.j )  d^  ...  dps  =  0(N), 


1 


DN(x1®...0AJ)dx1  ...  dXj  =  0(1) 


Thus, 


0  J 


cum (d-j  1(^1  )c*]2(^l )  *  •  •  * »  dj-| (*j)dj2(*j))  dA-| . .  .dx .  -  0(N), 


and  consequently 


cum(A.  ,  . . . ,  A.  )  =  N"J//2 


V  TJ 


1  rl 


♦  Ui)  ...  ♦  ,  (*,) 
0  1 1  7J  J 


cum  (d^  -j  ( )d^2(x^ ) ,  ...,  dj-j(Xj)dj2(Xj))dx^  ...  dx  j 
=  0(N"J/2+1), 

which  implies  the  asymptotic  normality,  q 


Suppose  f(x)  is  the  spectral  density  of  a  stationary  process  and 
{ffl(x)}  is  a  family  of  fitted  spectral  densities  which  are  parameterized 
by  e  e  0  c  Kr,  where  0  is  a  compact  set  in  IRr.  We  define  a  pseudo-true 
value  e  of  e  e  0  C  lRr,  by  a  value  which  minimizes 


D(f0,f)  - 


(log  det  f  (a)  +  tr  f(x)ftt(x)"  )dx 

0  u 


with  respect  to  e  e  0. 

ASSUMPTION  2.  The  fitted  model  fQ(x)  is  twice  continuously  differ- 

0 


entiable  with  respect  to  e  e  0. 


'It  *l». 


•4^*. 


», 


A 


ASSUMPTION  3.  If  §  f  e*,  then  fQ(x)  t  f0*(X)  on  a  set  of  positive 
Lebesgue  measure.  The  matrix 


Me)  =  ' 
1  ~ 


■1  s2  1 

3030^-  [log  det  fg(x)  +  tr  f^ (x)  'f(A)]dA, 


is  nonsingular  for  all  e  ee,  and  M^.  =  M^(e). 

The  first  statement  of  Assumption  3  is  an  identifiaility  condition. 
In  Section  4  some  nonidentifiable  examples  will  be  given.  Then  we  have 
the  following  theorem. 

THEOREM  1.  Let  { Y( t ) }  be  a  q-dimensional  dyadic  stationary  process 
with  mean  zero  and  the  spectral  density  f(x).  Suppose  that  Assumptions 
1-3  are  satisfied,  and  that  §  exists  uniquely  and  lies  in  Int  0.  Then 

(i)  lim  e  =  e  in  P, 

N-**> 


(ii)  the  distribution  of  the  vector  /N{P-§},  as  N  -*  <*>,  tends  to  the 

normal  distribution  with  mean  zero  and  covariance  matrix  M^1  V  M'1 ,  where 

V  =  (V.  }  is  a  r*r  matrix  such  that 
jm 


V.  =2 

jm 


f1  a 

tr[f(x)  jf- 

J0  aej 


(f9(x))-'f(x)  4-  tfe(x))_1]  !  dx 

-  m  -  - 


I  Hrr- 

a,b,c,d=l  -*0'  (  30j 


1  (3fJb,a)(X)  3fJd’c)(x) 


9=e  fabcd^’x,y^dAdlJ’ 


where  fib *a ^ ( X )  is  the  (b.a)-th  element  of  f  ( X ) ” 1 . 
0  “ 


Proof.  From  the  definitions  of  e  and  e,  we  have 


h  D<f9-IN>8=e  ■ 


v"-: 


Te  0(fe,f)e=5  =  o. 


Expanding  (2.14)  around  e,  we  have 


0  ■^D‘fvIN)  +ftf(s*)(9-i) 


where  0*  lies  on  the  straight  section  with  end  points  e  and  0,  and 

=  jheT  D(fe*’V' 


(2.15) 


(2.16) 


By  Lemma  1 ,  we  have 


&D(fg,IN)  -  0,  in  P 


Mf(§)  Mf(e)»  in  P  for  each  0  e  0- 

By  Assumptions  2  and  3,  absolute  values  of  eigenvalues  of  Mf(e)  have  a 
positive  lower  bound  for  all  9  e  0.  Hence  when  n  is  large  enough,  with  a 
probability  arbitrarily  near  to  one,  so  do  the  absolute  values  of  eigenvalues 


of  Mf(e).  By  (2.16)  we  have 


and  consequently 


§  -v  0,  in  P 


Mf(e*)  Mf(§),  in  P. 

Then  the  limiting  distribution  of  /N(e-e)  is  equivalent  to  that  of 

-M^4D(f5’IN) 

=  -M'1^  J1  ^  (log  det  fjj(A)  +  tr  f g (x ) ' 1 IN ( A ) } dx 
=  j1  £tr  fe(A)"1UN(A)  -  f ( A )) ] dx , 


(2.17) 


rm 


by  (2.15).  Again  applying  Lemma  1  to  (2.17),  we  have  completed  the 
proof,  p 

Remark.  If  the  true  Walsh  spectral  density  matrix  f(x)  =  fD(x), 

0 

the  pseudo-true  value  is  equal  to  the  true  value,  i.e,,  0=0  (see 

Hosoya  and  Taniguchi  (1982)). 


3.  MODEL  SELECTION  OF  WALSH  SPECTRAL  MODELS 


In  the  previous  section  we  assumed  that  the  order  of  the  unknown 

parameter  vector  of  the  Walsh  spectral  model  f _ ( X )  is  known.  However,  in 

the  actual  situation,  we  must  estimate  the  order  of  dim  9  =  r  from  the  data. 

Here  we  assume  that  the  process  { Y ( t ) >  has  the  Walsh  spectral  density  matrix 

f.  (a),  0  =  (e, ,  . . . , 6  ) 1 ,  where  0  is  an  unknown  parameter  vector.  (We 

~r  i  r  ~r 

use  suffix  r  to  stress  the  dimension.)  Then  we  fit  the  Walsh  spectral 

model  f .  (x),  0  <  k  <  L,  where  L  is  a  preassigned  upper  limit  to  the  order. 

§k  “  "  . 

We  determine  the  true  order  r  by  the  value  k  which  minimizes  the  following 

criterion: 


A(k)  .  k-y 


for  k  =  0,1 .... ,L, 


(3.1) 


where  C^  -*■  ®  and  C^/N  +  0  as  N  -*■  ®.  For  this  estimated  order  k  we  have 

THEOREM  2.  Suppose  that  all  the  assumtions  in  Section  2  for  f(A)  = 
f  (a)  and  f0U)  =  fQ  (a)  are  satisfied.  Then  lim^k  =  r,  in  p. 


Proof .  From  (?.16)  we  have 

/R(ek-ek)  =  '^f1(?k)/R  1  af-  [tr  f§  (^)_1Un(a)  -  f(A))]dA 

0  -  k  k 

which  tends  to  normal  by  Theorem  1.  Thus  we  have,  for  any  sequence  of 
positive  numbers  C^  -*•  ®, 

P[\\  >  CN]  =  o(l), 

where  ||  •  ||  is  the  Euclidian  norm.  Taking  CN  =  Vc^t  we  obtain 


(3.2) 


(3.31 


ft 


Expanding  around  8  =  e^,  and  noting  (3.4)  we  can  see  that 


3D(f0  *IN) 

=  D(fek,IN)  +  "  -k^ '  ^ 


V=  §k 


+  ^ik-ik^*  Rf{®k)(§k-?k> 


(3.5) 


Si  nee 


3D(f8k,IN) 
'  8®k . 


-V! 


,  =  $,  we  have 


o<fek -'n5  =  D<f|fc-!N>  -I(§\-Ik),¥-ek)(2k-§k)' 

As  first  step  we  show  that 

P(k  <  r)  -*■  0  as  N 


(3.6) 


(3.7) 


For  k  <  r,  we  evaluate 

/  \  ^  r—k ) Cm 

P1  =  P (a ( k )  <  A(r)j  =  P(D(fg  ,IN)  -  D(fg  ,IN)  <  - >. 

*  k  **  r 

Using  the  relation  (3.6),  the  above  probability  is  approximated  as 

p{D<f§k'V  -  D<fer>V 

(r-k)CN  ,  _  * 

<  N  +  2  (®k"fk),Mf(®k)(lk-®k) 

"  I(?r'  ?rJ'Mf(?r)(?r'  ?r)K  (3,8) 


Using  Lemma  1  the  left-hand  side  of  the  above  {•}  converges  to  D(fr  ,f)  - 

-k 


D(f-  ,f),  which  is  strictly  positive  for  k  <  r.  On  the  other  hand,  by  (3.4), 
~r 

the  right-hand  side  of  {•}  converges  to  zero  in  probability  which  implies  the 


•.  ■.  •- 


probability  0  as  N  -*■  ».  As  second  step  we  show 


P(k  >  r)  -*•  0  as  N  -*■  ». 


(3.9) 


We  have  for  k  >  r. 


(r-k)CN 

P2  =  P{ A( k )  <  A(r)}  =  P{ D(fg  ,IN)  -  D(fg  ,IN)  <  - -1} 

k  ^  r 


Using  the  relation  (3.6),  the  above  probability  is  approximated  as 


P(D(f-  ,iN)-D(f-..iN)-i(ek-?k)'V!k)(5k-!k> 

~  k  ~r 


+  ^(®r"  !r,,ftf^r^§r“  ?r}  < 


(r-k)C. 


(3.10) 


(3.11) 


Because  f-  (a)  -  f-  (x),  for  k  >_  r,  we  can  see  that 
~k  ~r 


D(f5  .iN)  -  D(f5  ,IN)  -  0. 
-  k  ~  r 


While,  by  (3.4),  we  can  see  that 

•  1  <?k  -  !k> '  V5k»?k  -  Sk>  +  1  -  5r>  'Sf  <?r>  (?r  ‘  ?rJ 


is  at  most  of  order  0p(^p-).  However  the  right-hand  side  of  {•}  in  (3.10)  is 
(r-k)CN 

- ^ — - ,  (r<k),  which  implies  P2  ->  0,  as  N  -*•  ».  Thus  we  have  completed 

the  proof.  □ 


VvVv.vVV.V.-.-.V.VV  V.  .  V.  .  A  A 
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4.  DETERMINATION  OF  THE  ORDER  OF  DYADIC  AUTOREGRESSIVE  MODELS 


* 

V 


3 


In  the  previous  sections  we  could  proceed  in  fairly  analogous  ways 
to  those  used  in  the  ordinary  stationary  processes.  However  if  we  consider 
finite  parametric  models,  for  examples,  dyadic  autoregressive  process  of 
finite  order  (DAR-process) ,  dyadic  moving  average  process  of  finite  order 
(DMA-process)  and  dyadic  autoregressive  moving  average  process  of  finite 
order  (DARMA-process) ,  then  there  exist  the  greatest  differences  between 
dyadic  stationary  processes  and  ordinary  stationary  ones.  That  is,  it  is 
known  that  these  DAR,  DMA  and  DARMA  are  equivalent  in  the  sense  that  DAR 
or  DARMA-process  of  finite  order  can  be  expressed  as  DMA-process  of  finite 
order  (see  Nagai  (1980)  or  Nanai  arid  Taniguchi  (1987)). 

In  this  section,  for  a  finite  order  dyadic  autoregressive  model,  we 
can  propose  a  simpler  order  determination  criterion.  Then  we  show  that  this 
criterion  gives  strong  consistent  order  estimate. 

A  q-dimensional  dyadic  stationary  process  (Y(t),  t  e  T}  is  called  a 
dyadic  autoregressive  process,  if  it  can  be  expressed  by 

![  A.Y(t©  j)  =  e  (t) ,  t  e  T,  (4.1) 

j=0  J 

where 

(i)  A.'s  are  qxq  matrices,  An  =  I  ,  and  p  =  2r  -  1,  where  r  is  a  non- 
negative  integer, 

(ii)  e(t),  t  e  T,  are  i.i.d.  random  vectors  such  that 


a 


3 


,V 

.0 


(iii) 


Ee( t)  =  0,  Ec(t)e(t) '  =  G  >  0, 


det  *(x)  f  0,  a.e.x , 


(4.2) 


(4.3) 


where  *  (x)  =  \  A.W( j  ,X) , 
j=0  J 


£ 

S 

vi 


If  (4.2)  and  (4.3)  hold  true,  then  the  Walsh  spectral  density  of  { V ( t ) }  is 


f(x)  =  <!>(  a)_1G{  <t>(A)-1}  * . 


(4.4) 


We  call  a  DAR  process  (4.1)  irreducible  if  there  does  pot  exist  such 


a  matrix 


2r-h 


which  satisfies 


*XU)  =  }  K,W(j,x), 


f(x)  =  ,  a.e.  X. 


(4.5) 


Especially,  for  an  irreducible  DAR  process  (4.1),  there  exists  a  tQ, 

2r_1  <  tn  <  2r  -  1 ,  such  that  A.  i  0. 
u 

For  an  irreducible  DAR  model  (4.1),  p  is  called  the  order  of  the  model, 
For  simplicity,  such  a  model  is  written  as  DAR(p).  Note  that  in  the  above 
definition,  the  order  of  the  model  (4.1)  is  defined  as  p  =  2r  .  l ,  not  as 
max{t:  f  0} .  The  advantage  of  such  a  definition  is  that  it  suits  to  the 

Walsh  spectrum  analysis,  and  is  convenient  for  estimating  the  parameters 
of  the  model.  To  see  this,  consider  the  following  two  scalar  irreducible 


DAR  model : 


X(t)  +  X(t®  1)  +  aX(t®  2)  =  e(t) , 


N 


m 


Y(t)  +  Y(t@l)  +  aY(  t  (T)  3)  =  e(t),  t  e  T, 

where  a  ^  0,  a  f  ±2,  e(t)'s  are  i.i.d.  with  Ee(t)  =  0,  Ee(t)^  =  It  is 

easily  seen  that  they  have  the  same  Walsh  spectral  density 

a^[l  +  W(l,x)  +  aW(2,x)] 

But  if  we  define  the  order  of  the  model  as  max{t:  A^O},  then  their  order 


may  be  2  and  3  respectively.  Obviously  such  a  definition  is  not  convenient 
for  Walsh  spectrum  analysis.  It  is  easy  to  see  that  these  two  models  are 


not  essentially  different.  For  a  q*q  matrix  A  =  (a..,  1  <j  ,  j  <_  q ) ,  denote 

|| A ||  =  .  ,  | a . . | .  To  determine  the  order  p  =  2r  -  1  of  the  irreducible 

•  * J  I  I  J 

model  (4.1),  we  suggest  the  following  criterion: 


2^  i  N  1 

'  L"(k)  X  "w  J0Vtt>v(t©  (2k  +  "))  • 


(4.6) 


where  Y(0) ,  . . . ,  Y ( N- 1 )  are  the  observations  of  the  model  (4.1),  N  =  2m  with 
m  positive  integer,  and  satisfies  the  following  conditions: 


1  im 

|\|-yco 


0 


and 


1  im 
N-k” 


log  log  N 


(4. 


Define 

rN  =  max{k^0:  l_N(k-l)>0,  l_N(k)<0}.  (4.8) 

where  LN(-1)  =  1  for  convenience.  We  can  use  r^  as  an  estimate  of  the  true 
value  r  of  the  model  (4.1).  We  have  the  following 


THEOREM  3.  If  the  model  (4.1)  is  irreducible  and  (i),  (ii)  and  (iii) 
are  satisfied,  then 

lim  r,,  =  r,  a.s„  (4.9) 


Proof.  Suppose  that  p  is  the  true  order  of  the  model  (4.1)  and 
p  =  2r  -  1.  According  to  Nagai  and  Taniguchi  (1987),  if  det{$(x)}  1  0, 
then  { Y ( t ) :  t  e  T}  is  a  DMA-process  written  by 

2r-1 

Y(t)  =  I  t  e  T.  (4 


10' 


Put  r(n)  =  EV(0)Y'(n).  By  (4.10)  and  the  condition  (ii),  it  is 
easily  seen  that  for  any  n, 

lljV  Y(t)Y'(t©n)  -  r(n)  ||  =  0(/log^g~) ,  a.s.  (4. 

as  N  ->■  <*>,  (e.g.,  Petrov  (1975)).  By  (4.10),  for  n  >  2r,  r(n)  =  0.  Thus, 


if  k  >  r,  then 


LN(k)  - 


as  N  -v  °°.  From  this  and  ^  «  ®,  it  follows  that  with  probability 


one  for  N  large. 


i-N(k)  <0,  k  >_  r. 


If  r  =  0,  the  theorem  is  proved. 


Now  assume  that  r  >  0.  We  have 


2r-1-l  ,  N-l 


Tim  l  II  jr  l  Y(t)Y(t©(2r_1  +n)) 
n=0  Nt=0  v  7 


2r_1-l 


l  ||r(2r-1+n)||2,  a.s. 
n=0 


We  proceed  to  prove  that 


Otherwise,  we  have 


2r_^_i 

l  ||r(2r'  +  n)  ||2  >  0. 
n=0 


2r-’-l 


f(x)  "  L  x  e  [0,1]. 


wrongs? 


20 


Put  h  =  2 r_  - 1 ,  x.  =  j/(h  +  l),  j  =  0,1,. ..,h.  We  know  that  for  all 
J 

l  <_h,  W U,x)  =  W(£,Aj)  for  x  e  [ Aj !  xj+-j )  *  From  this  it  1S  easily  seen 
that  f ( x)  takes  only  at  most  h  +  1  different  values,  say,  f (Aq)  , . . .  ,f )x^ 
By  (4.4),  G  >  0  and  <j>(x)  f  0,  it  is  easily  seen  that  f ( x  . )  >  0,  j  =  0,1, 

<J 


Hence  we  can  write 


G  =  G1/2(G1/2)’,  f(x.)  =  f'/c(x  )f'/c(x  )', 

vJ  J  J 


*  ,1/2 „  -1/2, 


j  -  0,1 ....  ,h.  Put 


'W(0,x0)  W(1 ,XQ)  ...,  W(h ,XQ) 
WlO.X^  Wd.X^  ....  W(h  ,X  1 ) 


W(0,xh)  W(l,xh)  ....  w(h,xh) 


Then  H^+1Hh+1  =  ( h+1 ) I h+1 .  Thus  the  matrix  equation 


<Vi®v  : 


G-1/2f'/2(x0) 


G-1/2f,/2(xh) 


has  a  unique  solution  (Bq,...,B^),  where  B^.'s  are  all  q*q  matrices.  Fn 
(4.17)  we  can  see  that 

h  h  x 

(  l  B  W(£,x,)  G(  l  B  W(i,x.)  1  =  f ( X • ) ,  j  =  0,1 . h,  (4 

'o=n  *  J  o=n  x  J  j 


which  implies 


n ( X ) Gn ( X ) 1  =  f ( X) ,  x  e  [0,1] 


n(A)  =  l  B  W(£,X). 


£=0 


By  Nagai  and  Taniguchi  (1987)  there  exists 

2  r" 1  - 1 

*i(a)  =  l  K  W(£,A) 
£=0 

such  that 


4>lU)n(x)  =  Iq,  a.e.  x. 


Thus  we  have 


f(x)  =  *1(x)“1G{$1(x)"1}' , 


a.e. 


A  5 


( 


which  contradicts  our  i rreducibi 1 i ty  assumption.  Now  (4.15)  has  been 
proved.  By  (4.14),  (4.15),  (4.6)  and  (4.7),  with  probability  one  for 
large  N, 


LN(r  -  1 )  >  0,  r  >  0. 

Noting  (4.13)  and  (4.23),  with  probability  one  for  large  N,  we  have 


rN  "  r* 


Remark „  The  following  scalar  process  { Y ( t ) ;  t  e  T)  is  a  reducible 
DAR  process : 

X(t)  +  X(t®l)  +  X ( t® 2 )  -  X(t®3)  =  e(t),  t  €  T, 

2  2 

where  e(t)'s  are  i.i.d.  with  Ee(t)  =  0  and  Ec(t)  =  a  . 

*(A)  =  1  +  W(  1  ,  A)  +  W(  2  ,  A)  -  W(  3  ,  A)  , 


Then 


5.  DETECTION  OF  SIGNALS  FOR  DYADIC  STATIONARY  PROCESSES 


In  this  section  we  shall  consider  a  signal  detection  model  for  dyadic 
process  of  finite  order,  and  show  that  this  model  is  equivalent  to  a  dyadic 
moving  average  model.  Then  we  can  apply  the  results  in  Sections  3  and  4  to 
our  model.  That  is,  we  can  determine  the  order  of  the  signal  dyadic  process 
from  the  data. 

Let  { Y ( t ) }  be  a  q-dimensional  dyadic  stationary  process  defined  by 

P 

Y(t)  =  l  A(j)S(t@j)  +  5(t),  (5.1) 

j=0 

where  S(t)  is  an  r-dimensional  dyadic  stationary  signal  process,  and  c(t)  is 
a  q-dimensional  dyadic  stationary  noise  process,  and  A(j)  are  qxr-matrices . 
Here  we  assume  that  { S ( t ) >  is  an  r-dimensional  DARMA(s ,h) -process  defined  by 

s  h 

l  B(j)S(j@t)  =  l  C(j)U(j0t),  t  =  0,1 .  (5.2) 

j=0  j=o 


where  l U( t ) }  is  an  r-dimensional  white  noise  process.  Also  U(t)}  is  a 
q-dimensional  DARMA(£ ,m)-process  defined  by 


l  D(jk(j0t)  =  l  F ( j ) V ( j 0 1 )  ,  t  =  0,1,...,  (5.3) 

j=0  j=0 

where  { V ( t ) }  is  a  q-dimensional  white  noise  process  which  is  independent  of 
{ U( t) } .  We  assume  that  all  the  coefficients  { A( j ) } ,  (B(j)},  { C( j ) } ,  { D( j ) } 
and  (F(j)}  are  completely  specified  by  an  unknown  parameter  vector  e  with 
dimension  k,  and  that 
s 

det{  I  B( j )W( j  ,x) }  f  0, 
j=0 

£ 

det{  \  D(j)W(j,a)}  f  0,  for  all  0  <_  X  _<  1 . 
j=0 


(5.4) 


Denote  by  Walsh  spectral  representations  of  { Y ( t ) } ,  { S ( t ) } ,  U(t)}, 

(U(t)}  and  tV(t)}  by  Y(t)  =  /Jw(t,x)dZy(x) ,  S(t)  =  /Jw(t,x)dZ$(x) , 

?(t)  =  /Jw(t,x)dZg(x) ,  U(t)  =  /Jw(t,x)dZu(x)  and  V(t)  =  /Jw(t,x)dZv(x) , 
respectively.  The  relation  (5.1)  can  be  written  as 

dZy(x)  =  TA(x)dZs(x)  +  dZ^(x) ,  (5.5) 

briefly,  where  ta(x)  =  jJ=QA( j )W( j ,x) .  While  we  can  have 

dzs(x)  =  4'B(x)”1H'c(x)dZ|J(x),  (5.6) 

dZ5(x)  =  H'D(x)"1M'F(x)dZv(x),  (5.6) 

where  fg(x)  =  ^=QB(j)W(j  ,x) ,  *c(x)  =  Xj=0c(j )W( j  ,x) ,  yq(x)  =  £*=QD( j )W( j  ,x) 
and  fF(x)  =  X j=0F( j )W( j ,x ) .  Finally,  we  get 

dZY(^)  =  ’*'A(^)TB(x)"14'c(x)dZu(x)  +  T[)(x)"1'}'F(x)dZv(x) 

i  i  rdz„(x)- 

=  0A(xhB(x)  Tc(x):  vD(x)  SF(X)]  .  (5.8) 


By  Nagai  and  Taniguchi  (1987),  under  (5.4),  we  can  get  the  following  repre¬ 
sentations 

1  s’ 

vB(x)  =  l  B( j )W( j ,x) 

j=0 


and 


i' 


Tn(x)'1  =  l  D( j )w( j ,x) , 

u  j=0 


where  s' =2  -  1,  -  2  -  1,  a  and  b  are  the  minimum  nonnegative  integers 
which  satisfy  s  <_  2a  -  1,  £  <  2^  -  1,  respectively.  Thus  it  is  easy  to  see 
that  4'a(x)4'b(x)"\(,(x)  and  ^(x)-1  ^(x)  can  be  written  as  finite  linear 


combinations  of  Walsh  functions 

dl  dl 

l  G(j)W(j,x),  l  H( j )W(  j ,x) , 

j=0  j=0 

respectively.  So  the  relation  (5.8)  implies  that  our  model  (5.1)  is  equiva 
lent  to  a  DMA-process  of  finite  order.  We  can  assume  that  our  { Y ( t) >  has 
the  finite  DMA  type  Walsh  spectral  density  matrix  fQ(x).  To  determine  the 
order  of  dim  0,  we  can  use  the  criterion 


A(k)  =  Dtf^.  IN)  + 


given  in  (3.1).  Of  course  we  can  use  A(k)  to  determine  the  number  of 
signals . 


6.  TEST  OF  HYPOTHESIS  FOR  LINEAR  RESTRICTION  OF  PARAMETERS 


Let  T Y ( t ) }  be  a  scalar-valued  dyadic  stationary  process  with  Walsh 
spectral  density  fQ(A)  depending  on  an  unknown  parameter  e  =  (e^,  . . . ,  e ^ . 
We  assume  that  {Y(t)>  satisfies  all  the  assumptions  in  Theorem  1.  The 
first  problem  is  to  test  a  composite  hypothesis  H^:  =  e2g?  against 

Hi  §2  ^  §20  *  §  -  (§i  »§2^ »  §  i  ”  ( ®  i »  **‘,0£^*  §2  _  £^1  *  • « • »  ®  p ) 

and  e^g  =  (e£+^  g»  ...,6  q),  a  specified  vector  and  (Si*?20>  e  Int  0. 

Although  we  do  not  assume  tne  Gaussianity  of  (Y(t)},  we  can  formally  make 
the  following  log-likelihood  ratio  criterion 


ZlogL  =  N{D(f(g^),  IN 


}  "  D(f(|r02O)’lN)>’ 


(6.1 


where  e1  =  (ej,^)  is  the  quasi-maximum  likelihood  estimator  for  e  under  H, 
and  e-j  is  that  for  under  HQ.  Put  v  =  ;ft(0-0),  w  =  vTT(0-|-0-|)  and  u1  = 
(w1, (J) 1 ) .  Expanding  in  a  Taylor  expansion  around  e,  we  have 


-G  = 


(u  -  v)  (l  +Op(D) 


1  \ ,  a  D^f§’  V 


2"  (u  -  v) 1  9639 


1, 

2  ^ 1  9090 


(u-v)(l+op(l)J 

=  |-(u  -  v)’Mf(u  -  v)  *  (l  +  op(  1 ) )  • 


(6.2 


From  Theorem  1  we  have 

v  (6'3 

Similarly  we  have 


u*  V^VV^V1))- 


(6.4 


Jr1  *4* 


_ 


where 


!il  0 
0  0 


.U  =  In  +  0(N~  ) , 


(®1  •§20)*  N  "  11 


From  (6.2),  (6.3)  and  (6.4)  we  have 


,  3D(fe,I  )  ,  .  3D(fe,Iw)  , 

_G  =  _/N  [Mf  -  Lf]Mf[Mf  -  Lf]/iT  ^  (1+opd)) 


1  _  ^(  »  In)  -1 


3D(fg  »Im)  . 

— sS— “-(l+OpCll) 


(6.5) 


Here  we  put  the  following  assumptions: 


Assumption  4.  The  process  { Y ( t ) )  is  a  scalar  linear  dyadic  stationary 


process  represented  as 


Y(t)  =  l  A.e(t®  j), 
j=0  J 


(6.6) 


where  V._Q|A.|  <  °°  and  the  e(t)'s  are  independent  random  variables. 

J  u  J 

Assumption  5.  The  unknown  parameter  0  of  f  (x)  is  innovation-free,  i.e. 

“*  0 


4  l!lVxr'fe„(x),dx  = 


(6.7) 


(See  Hosoya  and  Taniguchi  (1982).) 


LEMMA  2.  Suppose  that  Assumptions  1-5  are  satisfied,  and  that 
Ij=0AjW(j,x)  f  0  for  all  X  e  [0,1].  For  an  innovation-free  parameter  g  we 


'N  9D^  W  £ 


N(*p.Mf). 


(6.8) 


Proof.  Using  a  similar  argument  to  Hosoya  and  Taniguchi  (1982),  we 


can  see  that 


If  3ffl (X )  f  g  (  X  ) 


0i  36r 


f^U.X.yJdAdy  =  0, 


for  j,m  =  where  f^(*)  is  the  fourth  order  cumulant  spectral  density 

Putting  <|>.(x)  =  r|—  fQ(x)  in  Lemma  1,  we  have  the  desired  result,  n 
j  99j  e 


Applying  Lemma  2  to  (6.5)  we  have 

THEOREM  4.  Suppose  that  Assumptions  1-5  are  satisfied.  Then  the  dis- 
tribution  of  -G  under  Hq  tends  to  x  (p-  0  as  N  +  ®. 

Now  we  consider  a  more  general  test  of  hypothesis. 

Hq :  Be  =  u2Q  against  H:  Be  f  u2Q, 

where  B  is  a  (p  -  £)*p  matrix  with  rank  B  =  p  -  and  u2Q  =  (u^+1  Q,  . . . ,  up  Q 
Then  there  exists  an  £xp  matrix  A  such  that 


(  B  )-  *  [ll)- 


where  det(g)  f  0.  Let  e  be  the  quasi-maximum  likelihood  estimator  of  e  e  0, 
then  u(e)  =  0.  Then  the  likelihood  ratio  criterion  of  testing 


Hq:  u2  =  u2q  against  H:  u7  f  u 


is  given  by 


G  =  N{D(f 


(G1  ,G2)  •  IN)  "  D(f(Gru20)’  JN 


2  '  20 


.  IJ>. 


(6.9) 


where  u^  is  the  quasi-maximum  likelihood  estimator  of  u^  under  Hq.  Then 


we  have 


THEOREM  5.  Suppose  that  Assumptions  1-5  are  satisfied.  Then  the  dis- 

2 

trihutions  of  -G  under  Hn  tends  to  y  (p-  &)  as  N  ->  ». 
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